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ABSTRACT 

I consider the physics of gravitational instabilities in the presence of dynamically important radiation pres- 
sure and gray radiative diffusion, governed by a constant opacity k. For any non-zero radiation diffusion rate 
on an optically-thick scale k , the medium is unstable unless the classical gas-only isothermal Jeans criterion 
is satisfied. If the radiation acoustic sound crossing timescale on a scale k~ l (t r ) is less than the geometric 
mean of the dynamical and radiative diffusion timescales (tats and ?dy n ), then diffusion is "slow." In this limit, 
although the dynamical Jeans instability is stabilized by radiation pressure on scales smaller than the adiabatic 
Jeans length, on these same spatial scales the medium is unstable to a diffusive mode. In this regime, neglect- 
ing gas pressure, the characteristic timescale for growth of this mode is independent of spatial scale and given 
by (3kCj)/(47tGc), where c s is the adiabatic sound speed. This characteristic timescale is that required for a 
fluid parcel to radiate away its thermal energy content at the Eddington limit, the Kelvin-Helmholz timescale 
for a radiation pressure-supported self-gravitating object. In the limit of "rapid" diffusion — defined by the 
inequality t r > (t&g fdyn) I//2 — radiation does nothing to suppress the Jeans instability and the medium is dy- 
namically unstable unless the gas-only Jeans criterion is satisfied. I connect with treatments of Silk damping in 
the early universe. I briefly discuss several astrophysical applications, including photons diffusing in regions of 
extreme star formation (starburst galaxies & pc-scale AGN disks), and the diffusion of cosmic rays in normal 
galaxies and galaxy clusters. The former (particularly, starbursts) are "rapidly" diffusing and thus cannot be 
supported against dynamical instability in the linear regime by radiation pressure alone. The latter are more 
nearly "slowly" diffusing. I speculate that the turbulence in starbursts may be driven by the dynamical cou- 
pling between the radiation field and the self-gravitating gas, perhaps mediated by magnetic fields, and that 
this diffusive instability operates in individual massive stars. Appendices contain a more detailed treatment of 
radiation transport and consideration of uniform rotation in the background medium. 

Subject headings: Instabilities, Hydrodynamics, Radiative Transfer, ISM: Cosmic Rays, Stars: Oscillations, 
Galaxies: Starburst 



1. INTRODUCTION 

Standard treatments of the Jeans instability assume the 
medium is homogeneous and isotropic and governed by a 
barotropic equation of state. Employing the "Jeans Swindle" 
so that the Poisson equation is satisfied in an ad hoc way with 
no background gradients in density, the dispersion relation 



= c z „k l -4-nGp 



(1) 



follows from a linear analysis. Here, c g is the gas sound 
speed and p is the mass density. The Jeans instability is long- 
wavelength; for scales larger than the Jeans length, 



2nk, 



:Aj= C ,(VGp) 1/2 , 



(2) 



the system is dynamically unstable under the action of a per- 
turbation to the density and the attending increase in the gravi- 
tational potential (Jeans 1902, 1928; e.g., Binney & Tremaine 
1987). On scales smaller than Aj the medium responds to a 
compression with a restoring pressure force. Equation (fJJ can 
be obtained by equating the acoustic sound crossing timescale 
on a scale Aj with the dynamical timescale. Equivalently, 
equation (O may be read as expressing the fact that for stabil- 
ity, the total thermal energy of the medium within a volume Aj 
must exceed the gravitational potential energy. When the sys- 
tem is unstable, it continuously and spontaneously transitions 



to states of lower total energy by liberating thermal energy 
(e.g., Chandrasekhar 1961). 

The purpose of this paper is to understand how the classi- 
cal gas Jeans criterion is modified by radiation and to ask in 
which astrophysical environments such a modification might 
be important. Although there are many treatments in the lit- 
erature of both the Jeans instability (e.g., Jeans 1928; Ledoux 
1951; Chandrasekhar 1954, 1958; Mestel 1965; Lynden-Bell 
1966) and the physics of radiating flows (e.g., Mihalas & Mi- 
halas 1984 and references therein; Spiegel 1957; Kaneko et 
al. 1976; Bisnovatyi-Kogan & Blinnikov 1978, 1979; Mihalas 
& Mihalas 1983; Dzhalilov et al. 1992; Zhugzhda et al. 1993; 
Arons 1992; Bogdan et al. 1996; Gammie 1998; Kaneko et 
al. 2000; Blaes & Socrates 2001, 2003; Socrates et al. 2005), 
there has been relatively little work on self-gravitating envi- 
ronments where radiation might play an important dynamical 
role (however, see Kaneko & Morita 2006, Vranjes & Cadez 
1990; Vranjes 1990). 

Perhaps the first and most familiar treatment of self- 
gravitating radiation pressure dominated media was carried 
out by Silk (1967), (1968), and then extended by Peebles & 
Yu (1970) and Weinberg (1971), in the context of acoustic 
wave damping of primeval fluctuations by radiative diffusion 
— "Silk Damping" (see also Hu & Sugiyama 1996; Dodelson 
2003). However, these works focus specifically on the damp- 



ing rate of acoustic fluctuations and the generation of entropy 
and did not delineate how the Jeans criterion is modified on 
scales larger than the gas-only Jeans length when radiation is 
dynamically dominant and diffusing. They also do not discuss 
the physics of slow non-dynamical diffusive modes. 

This paper is motivated by astrophysical systems where 
self-gravity and radiation are essential. These include sites 
of extreme massive star formation such as compact star- 
burst galaxies and the parsec-scale disks or obscuring "torii" 
thought to attend the process of fueling active galactic nu- 
clei. These environments are marked by high radiation energy 
density and high gas density, as well as optical depths to their 
own dust-reprocessed infrared radiation that may significantly 
exceed unity ($3] and, e.g., Pier & Krolik 1992; Goodman 
2003; Sirko & Goodman 2003; Thompson et al. 2005 [TQM]; 
Chang et al. 2006). In each of these systems radiation pressure 
can be comparable to gravity and the associated photon en- 
ergy density is rivaled only by the energy density in turbulence 
and, potentially, the contributions from cosmic rays and mag- 
netic fields. ' The very high radiation energy densities in these 
systems led TQM to propose a theory of marginally Toomre- 
stable radiation pressure supported starburst and AGN disks 
(see also Scoville et al. 2001 & Scoville 2003). Additionally, 
this analysis may be of some interest for the stability of in- 
dividual massive stars and for self-gravitating media whose 
pressure is dominated by cosmic rays, either in normal star- 
forming galaxies or galaxy clusters. 

In §|2j I present a simple linear analysis of the gravita- 
tional instability. Appendix lAl discusses uniform rotation in 
the background medium. More detailed treatments of radia- 
tion transport are considered in AppendixlBl(see also Kaneko 
& Morita 2006). In §[3] I discuss the relevance of the results 
derived for a number of astrophysical environments and Sj4] 
provides a summary. 

2. GRAVITATIONAL INSTABILITY WITH RADIATION 

Here, I describe the simplest non-trivial treatment of the 
Jeans problem with radiation pressure and diffusion that cap- 
tures the physics needed for a qualitative understanding (see 
AppendixlBlfor a more detailed treatment). The equations ex- 
press continuity, momentum and total energy conservation, 
self-gravity, and equilibrium optically-thick radiative diffu- 
sion (e.g., Mihalas & Mihalas 1984). They are 



| + V.(pv) = 0, 

d\ 1 

_ + v-Vv = — VP-V$. 

at p 



and 



dU 

— + v-W + (t/+F)V-v = 
at 

V 2 $ = 47rG,o, 



c 
-- — Vu r . 
inp 



-V-F, 



(3) 

(4) 

(5) 
(6) 

(7) 



Here, U = u g + u r and P = p g + p r are the total internal energy 
density and pressure, and the subscripts r and g refer to the ra- 
diation and the gas, respectively. The radiation pressure force 

1 For a recent assessment of the strength of magnetic fields in starburst 
galaxies, see Thompson et al. (2006). For a discussion of cosmic ray feedback 
in galaxies, see Socrates et al. (2006). 



kF/c is contained in the VP term in equation (|4j. In addi- 
tion, p g = pk B T /m p , u r = aT 4 = 3p r , u g = Pg/i'f— 1), p is the 
gas mass density, 7 is the adiabatic index of the gas, F is the 
radiative flux, and k is the opacity. 

For simplicity, I take k constant and I do not distinguish be- 
tween the Planck-, flux-, and Rosseland-mean opacities. The 
above equations also neglect the time-dependence of the ra- 
diation field and they assume that the radiation and gas tem- 
peratures are exactly equal (see Appendix ITU see also Miha- 
las & Mihalas 1984; Gammie 1998; Blaes & Socrates 2003; 
Kaneko & Morita 2006). Because the Eddington approxima- 
tion has been made, the effects of photon viscosity have been 
neglected (e.g., Weinberg 1971; Agol & Krolik 1998). In ad- 
dition, in considering extreme star formation environments 
where radiation is reprocessed by dust, the above equations 
neglect the two-fluid nature of the coupled dust-gas system; 
that is, they assume perfect collisional and energetic coupling 
between the dust and gas (see §0. Finally, no terms repre- 
senting sources of optically-thin radiative heating or cooling 
are included. 

I consider perturbations of the form q — > q + 5qexp[ik • x- 
iu)t], keep only linear terms, employ the Jeans Swindle, 
and take the medium and radiation field as homogeneous, 
isotropic, and in radiative equilibrium: p = const, P = const, 
U = const, v = V-F = F = 0. The perturbed equations are 

-iojSp+ipk-S\ = 0, (8) 

-iujS\ + ik(SP/p) + ikS^ = 0, (9) 

-iujSU + (U + P)ik-8\ + ik-8F = 0, (10) 

-k 2 5<S>-4-TrG8p = 0, (11) 

ikSu r + (3K P /c)SF = 0. (12) 



and 



The thermodynamic perturbations to the total pressure and en- 
ergy density are 



op 
dp 



. dP 

8p+ of 

. dU 

, 6p+ w 



ST 



5T 



(13) 



(14) 



Note that (dp r /dp)\ T = (du r /dp)\ T = so that (dP/dp)\ T = 
(dp g /dp)\r = c\ — that is, only the gas makes a contribu- 
tion to the total isothermal sound speed, (dP jdp\j) l l 2 = cj, in 
equation ( fT3l . The perturbation to the radiation energy den- 
sity is written as 



Su r = 4aT 3 5T=A5T, 



(15) 



where the last equality defines A = du r /dT. Combining these 
thermodynamic relations with the perturbation equations, one 
finds that 



0=(uj 2 -c 2 T k 2 + 4irGp) ( 1 + — )+(c 2 T -c 2 )k 2 , (16) 
where 

rk 2 / A \ 

(17) 



. _ ck 2 f A 
3np \Cy 

is the diffusion rate on a scale k , Cy = (dU/dT)\ p is the 
total specific heat, c 2 = (dP/dp)\ s is the square of the adiabatic 



sound speed for the gas and radiation, S is the total entropy, 
and the identity 



OP 
dp 



dP 

J~p 



dP 
dU 



U + P 

P 



dU_ 

~dp~ 



(18) 



has been employed. Expanding the dispersion relation and 
combining terms, equation ( fT6] l becomes 

u? + iuj 2 £b-Lo [c 2 s k 2 -4-TrGp] -iu> [c 2 T k 2 -4-TrGp] =0. (19) 
Note that the dimensionless ratio 

A _ du r ( dU 
C^"dTydf p J \' ' 4u r/ 

that appears in equation ( fTTI i for ui approaches unity in the 
limit u s /(4u r ) — ► 0, and zero in the limit u g /(4u r ) — ► oo. 
Therefore, as u r /u g — + 0, c 2 
equation ( [T91 reduces to 



(20) 



c; =74 and (A/C v ) -> 0, 



uj[jc z T k z -4-TrGp] = 



(21) 



in the gas-pressure-dominated limit, fully analogous to the 
classical Jeans criterion in equation (Q}, but includes the en- 
tropy mode lu = (e.g., Lithwick & Goldreich 2001) and ex- 
plicitly contains the adiabatic gas sound speed c 2 = ^c\. In 
the opposite, radiation-pressure-dominated limit, u r /u g — > oo, 
c 2 — > c^ r = (4p r /3p) ^> c\. Neglecting gravity, the disper- 
sion relation for radiation pressure acoustic waves under the 
assumption of optically-thick equilibrium radiative diffusion, 
u)(oj 2 + iujui - c 2 r k 2 ) « 0, is obtained from equation ( 119) . 

2.1. Dimensionless Numbers 

Three dimensionless numbers determine the character of 
the modes admitted by equation (fT9l . The first measures the 
importance of gas pressure alone in supporting the medium 
on a scale k , 

(22) 



& T =c$.k i /(4nGp). 



The criterion -dj > 1 is the classical gas-only Jeans criterion 
for gravitational stability (cf. eq. J2j); $7- is the "isothermal 
Jeans number." The isothermal Jeans length follows by taking 

Aj. r = 2n/h, T = 2ttc t /{AttG P ) 112 . (23) 

The second and third dimensionless ratios combine to deter- 
mine the importance of radiation pressure. The first is 

$ s = c 2 k 2 /(4TTGp), (24) 

the "adiabatic Jeans number," analogous to $7-, but which in- 
cludes the contribution from radiation pressure. Taking i? s = 1 
defines the adiabatic Jeans length: 

Aj, s = 271-/**,, = 2ttc s /(4ttGp)^ 2 . (25) 

The second ratio is 

_ ck 2 f A\ 1 

a measure of the diffusion rate. The limits of rapid (x>l) 
and slow (x <C 1) diffusion are considered in Sections l2~2l and 
12.31 respectively. 

In analogy with the classical Jeans criterion, one might 
guess that if d s is larger than unity, then in the limit of slow 
diffusion the medium is stable. This turns out to be false, as 
I show in 32.31 In fact, if dj < 1 on a scale k , then the 
medium is unstable regardless of i} s . 

Using £ 2 = oj 2 /{4ttGp), and the definitions for dj, $,, and 
X, equation $1% can be written as 

£ 3 + ix£ 2 -«0*-D-&(0r-l) = O. (27) 



2.2. Rapid Diffusion 

In the limit of rapid diffusion (\ S> $. s , $r, 1), the three roots 
of equation d27j l are 



2x 



and 



^w-r'X + i 






(28) 
(29) 



to first order in \ x . When #r > 1, the roots in equation 
(f28b correspond to stable radiation- and gravity-modified gas 
acoustic waves. For large dj, these modes propagate at the 
isothermal sound speed of the gas; large \ ensures isother- 
mality. In the limit dj — > and \ 3> $s> equation d28l is sim- 
ply £ i=s ±/ and the medium is dynamically unstable. This is 
the classical gas-only isothermal Jeans instability. Note that 
the limit of rapid diffusion in equation d28| > is distinct from 
the high-A limit, because at high-A gravity, which dictates sta- 
bility/instability, disappears. However, to make an apposite 
comparison with the literature I take the high-A limit and for 
the acoustic modes I find that 

3/V " 



2k 
±c T k—i — u r 
3c 



1 + - 



Au r 



(30) 



in agreement with Blaes & Socrates (2003) (their eq. [62]). 
Equation d29l corresponds to the purely damped radiation dif- 
fusion wave. 

2.3. Slow Diffusion 
In the limit of slow diffusion (x <C #,, #r, 1), 



and 



£ « ±(0, 



£' 



2 



#,-1 



-'X 



•d T -l 
tf. s -l 



(31) 



(32) 



If $ s > 1, equation OTb corresponds to two stable damped 
gravity-modified radiation acoustic waves. 2 For # s 3> 1 and 
$s ^S> $t, the damping rate for these radiation acoustic waves 
is simply xj2. Conversely, when d s < 1 (and, thus, dj < 
1) — that is, on scales larger than the adiabatic Jeans length 
(eq. l25ll ) — the medium is dynamically unstable to the Jeans 
instability: £ — ► ±j, 

Equation d32b is key. It says that there is an intermedi- 
ate range in spatial scale k' 1 , larger than isothermal Jeans 
length (eq. |23l ) and smaller than the adiabatic Jeans length 
(eq. l25l ). that is always unstable. For negligible gas pressure 
($r — * 0), it is precisely when the adiabatic Jeans number #, is 
greater than unity and the dynamical Jeans instability is sup- 
pressed in equation ( TJTl i that the diffusive mode in equation 
( l32t is unstable. Even for arbitrarily large d s and small x, if 
the classical gas-only Jeans criterion indicates instability — 
that is, if -&t < 1 — then the medium is unstable. 

In a highly radiation pressure dominated medium with 
i9t "C 1 <C $ s , this diffusive mode grows at a rate 

MG ° (33) 



10 ■ 



3 k 



2 Note that the real part of £ is modified by \ at the level £ m ±($.j — 
1 ) ' /2 =F (X 2 /2) Wt - <tfr )(A ~ 1 r 3/2 in equation {31} if the second-order term 
in x is kept. 



in the high-£ limit, independent of spatial scale. This ex- 
pression is easy to understand as the rate at which a self- 
gravitating fluid parcel radiates its total thermal energy con- 
tent (e ~ (4/3)m,) at the Eddington limit (e = AirGpc/ k): it 
is the inverse of the Kelvin-Helmholz timescale: t^ ~ e'/e ~ 
4ttGc/(3kc2,). 

Equivalently, the only terms from the Euler and energy 
equations that contribute to this branch of the dispersion 
relation are the approximate equalities kSP/p s=s — k<5<E» and 
(U+P)k-S\ ps -k • 5F (cf. eqs.[9]&[T0]). Combining the energy 
equation with the continuity equation, {5p/p) = (k-S\)/u), and 
assuming that the medium is radiation pressure dominated, 

>c*\ rspA (i_p\ 

. K P/ \ 5 P J \ 4 PrJ 

On the other hand, the approximate equality kdP/p w -k(5$ 
implies that 

Combining equations 041 and (1331 ), one finds precisely equa- 
tion d33l . Substituting into the continuity equation, I find that 

k .to =i *p. = i(*^)(£), (36) 

up 4p r \ upc J \4-p r J 

which relates the density and velocity perturbations. 

Therefore, in a ra diatio n-pressure dominated medium with 
slow diffusion (see j$!2.51 l the characteristic time for collapse 
on scales smaller than the adiabatic Jeans length is ?kh, not 
the dynamical timescale. Although the growth timescale be- 
comes long as k and c s , r become large, for \ ¥ the medium 
is never formally stable if dj < 1 . Additionally, the existence 
of this instability does not require a pure radiation-only gas 
with adiabatic index of 4/3. Written another way, the growth 
timescale for this diffusive instability at large k is 

?KH ~ ?diff (fdyn Ar) (37) 

in an optically-thick, slowly-diffusing, radiation pressure 
dominated medium, where t r = (c sr k)~ l is the radiation pres- 
sure acoustic sound crossing timescale, and t&g = 3np/ck 2 and 
fdyn = (47rG/9) -1 / 2 are the diffusion and dynamical timescales, 
respectively. 

2.4. Criterion for Existence of the Radiation Acoustic Mode 

If the adiabatic Jeans number is larger than unity (i?, > 1), 
then when the radiation acoustic sound crossing timescale on 
a scale k~ l is shorter than the diffusion timescale across that 
same spatial scale, the radiation acoustic mode can be sup- 
ported by the medium. Thus, there is a critical diffusion rate 

defined by the rough inequality Xc\ > i?, , for which the ra- 
diation acoustic mode exists. This criterion on the diffusion 
rate \ can be obtained from an approximate solution to equa- 
tion (l27l in the limit $ s 3> 1 and $t — > 0, obtained by setting 
£ « for the radiation acoustic mode. I find that 
Xc,r _ 27T / c \ 1 
6 



2.5. Criterion Defining "Rapid" & "Slow" Diffusion 



1 



,1/2 



(38) 



2<&J ° \ c s,rJ T k 

where rj = K,p(2irk~ l ) is the optical depth on a spatial scale 
2irk~ l . Thus, if x is small with respect to Xc.r, then the dif- 
fusion rate on a scale k~ l is small compared to c s , r k and the 
radiation acoustic mode can be supported. Conversely, for 
X ^ Xc.r, such a mode does not exist. Note that the criti- 
cal value Xc,r is decreased by non-zero $r (see eq. J3T1 ) and 
modified if $ s is larger than, but near, unity. 



Sections 12.21 and 12.31 distinguish between the limits of 
"rapid" and "slow" diffusion. The criterion that separates 
these two limits defines a critical diffusion rate x c , diff that can 
be estimated by setting the growth timescale for the unstable 
diffusion mode in equation (1321 equal to unity, the inverse of 
the dynamical timescale. When d s ^> 1 ^> $7-, 

Xc,diff«<V (39) 

For x > Xc, diff the medium is "rapidly" diffusing and for x < 
Xc, diff it is "slowly" diffusing. Alternatively, equation 
may be written as (cf. eq. [37]) 



1 



ck 2 (AnGp) 



1/2 



3 up 



r- 



t r : 



(fdiff f dyn) ' '■ 



(40) 



if the radiation pressure acoustic sound crossing timescale on 
a scale k~ l is less than the geometric mean between the dif- 
fusion timescale on that same spatial scale and the dynamical 
timescale, then diffusion is "slow" and the stability proper- 
ties of the medium are best described by JJ2.3I Conversely, if 
tr > (fdifffdyn) 1//2 , diffusion is "rapid" Q2.2b . This criterion is 
valid only at high-A: and in that regime is independent of spa- 
tial scale. Equations ( [381 and ( [391 imply that Xc, diff can be 
greater than Xc. r and therefore that even though diffusion is 
"slow," the radiation acoustic mode is not supported. 

2.6. Solutions to the Dispersion Relation 

The limits of fast and slow diffusion, the criterion separat- 
ing them, and the range of existence of the acoustic modes 
and their damping rates are illustrated in the solution to equa- 
tion d27b presented in Figures[T]and[2j which show the modes 
£ obtained for a wide range of x, at fixed -d s and dj- Increas- 
ing x while keeping § s and dj constant can be thought of as 
a continuous decrease in the opacity k at fixed k~ l . Open and 
filled circles show the real and imaginary part of £, respec- 
tively. Individual pieces of the various roots are labeled for 
comparison with equations (f28l> — (f32t. 

The left panel of Figure Q] shows a case with $ s > 1 and 
t?r < 1. The unstable mode is the only positive imaginary 
root. For large \ it is the dynamical Jeans instability: £ w 
±/(l -^r) 1 / 2 , whereas for small x it is the diffusive mode of 
equation d32l . The dotted line shows the approximation to 
Xc.r (eq. 1381 ). Because for the parameters chosen, Xc>- ~ 

'l/2 

2zV ~ #s, the dotted line denoting \cr also roughly corre- 
sponds to Xc, diff (eq. l39l ). For x < Xc,r> the radiation acous- 
tic modes are evident and modestly damped. For any x/0, 
the medium is unstable because dj < 1 . Note that the purely 
damped mode £ sa — j'x is off-scale for large x (eq. l29l ). Con- 
trast the left panel of Figure Q] with the right panel, which 
shows the same calculation, but with $r = 3/2 > 1. Because 
#7- > 1, the medium is stable for any x- For the parameters 
chosen, the adiabatic acoustic mode, which contains contribu- 
tions from radiation and gas, joins smoothly into the isother- 
mal gas acoustic mode at modest x- 

Figure [2] presents a similar calculation, but more radia- 
tion pressure dominated, with d s = 20 and -dj = 1/10 (left 
panel) and dj = 3/2 (right panel). The left panel of Fig- 
ure Q] is qualitatively identical to the left panel of Figure [2j 
but in the latter there is a clear separation between Xo (dot- 
ted line) and Xc diff (dashed line). The right panel of Figure 



f «(l> s - l)'/2 
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FIG. 1 . — Le/i Panel: Solution to equation (27) for "d s = 3 and $r = 1/2, for 10 < \ 5: 10 at fixed £ . Open and filled circles show the real and imaginary 
parts of the three roots £, respectively. Positive complex roots indicate instability. Because '&j < 1 on the scale chosen, the me dium is unstable any \ =/0. For 
large \ the instability is dynamical, whereas for small x the medium is unstable to the diffusive instability given by equation )32t . The dotted line shows the 
approximate solution for Xcr (eq. 1381 ). below which the gravity-modified radiation acoustic mode exists with £ £s ±(i? s — l) 1 ' 2 . For the particular parameters 
chosen Xc.r ~ Xc.diff- Right Panel: Same as the left panel, but for -dj = 3/2. Because i9j- is larger than unity, the gravitational instability is stabilized for all \. 
Note the transition from adiabatic (radiation plus gas) to isothermal (gas only) gravity-modified acoustic waves. 
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FIG. 2. — Left Panel: Same as Figure[T] but for D s = 20 and -&t = 1/10, for 1 < \ < 10 2 at fixed IT 1 . Again, because i9j- < 1 on the scale chosen, the 
medium is unstable for any v =/ 0. For X > Xc, diff ^ #s (dashed line) the instability is dynamical, whereas for small X < Xc. diff the growth rate for instability 
is x(l — $t)/($s — 1) (eq. 1321 ). The dotted line denotes Xc, r (eq. |38| ). Right Panel: Same as the left panel, but $j = 3/2. Because "dj is larger than unity, the 
medium is stable for all \. The dashed line denotes the critical \ above which the isothermal gas acoustic wave exists with £ » ±(i9r — l) 1 ' 2 . 



|2] again shows a case with $7- > 1, so that the Jeans insta- 
bility is stabilized for any \< but shows the separation in \ 
between the adiabatic acoustic mode, which is here highly 
radiation pressure dominated (x < XcAt,r), and the damped 
gravity-modified isothermal gas acoustic mode that exists for 



X 



> 



(d s -d T )/2{d T -\) l l 2 (dashed line). 



Although I do not plot it here, equations (|281>-(l32j> and the 
general solution to equation (1271 1. show that in the limit ^7- — > 
and d s — » the medium is dynamically unstable for any \- 



2.7. The Connection to Treatments of Silk Damping 

The calculation of the dispersion relation and the slow 
diffusion limit of JJ2.3I is most closely related to the non- 
relativistic calculation of Silk (1967) in the context of the 
early universe. His expression for the damping rate of adia- 
batic radiation-dominated acoustic modes is essentially equiv- 
alent to the damping term in equation dXJV without the "1" in 
the denominator and with -dj — > 0. The correct relativistic 
expression for the damping rate was subsequently obtained 
by Weinberg (1971). The qualitative difference here with re- 
spect to Silk (1967) and other versions of the derivation of 
the damping of acoustic modes in the early universe by ra- 
diative diffusion (see, e.g., Hu & Sugiyama 1996) is equa- 
tion d32] i. which shows that the medium is unstable to a slow 
diffusive mode precisely in the regime {d s > 1) where the 
medium is stable to the adiabatic Jeans criterion — that is, 
on scales below the adiabatic Jeans length (eq. |251 ). A cal- 
culation of the dispersion relation in the cosmological context 
(in an expanding background) and an evaluation of the im- 
portance of this mode is in preparation. If it grows, it should 
do so very slowly on a timescale many times the dynamical 
timescale, and it should be purely non-adiabatic. Although I 
have not done a relativistic calculation here, taking c 2 ,. = c 2 /3 
in equation (l33l l. I find that co jw /47tG/(kc). This is the inverse 
of the Kelvin-Helmholtz timescale for a relativistic radiation 
pressure supported self-gravitating object. 3 The characteristic 
timescale for growth is of order ?kh ~ 10 16 (k/0.4 cm 2 g _1 ) s, 
which is of order thousands of times longer than the age of 
the universe at decoupling. Because the total matter density 
(which sets the gravitational driving term, the numerator of 
eq. 1331 ) is roughly ten times the baryon density (which sets 
the scattering timescale, the denominator of eq. [33)), one ex- 
pects a more careful derivation to yield a timescale shorter by 
this ratio. 



2.8. Extensions 

Motivated by Chandrasekhar (1961), Appendix lAl contains 
an analysis analogous to §|2] but including uniform rotation in 
the background medium. As in the case with only gas pres- 
sure, modes propagating at right angles to the angular mo- 
mentum vector are stabilized by rotation if the angular veloc- 
ity (51) is large enough that ft 2 > irGp (see also Goldreich & 
Lynden-Bell 1965). 

Appendix iBl accounts for the time-dependence of the radi- 
ation field and the possibility of imperfect energetic coupling 
between the radiation field and the gas. For the parameters 
appropriate to the astrophysical applications discussed in §|3] 

3 Note the correspondence with the Salpeter timescale for black hole 
growth. 



these factors are largely unimportant for the qualitative stabil- 
ity properties of the medium. This follows from the fact that 
the characteristic frequency for energetic coupling betw een 
the radiation and the gas, u>± sw Kpc(u r /u g ) (see eq. IIB8I and 
surrounding discussion, as well as Bogdan et al. 1996; Blaes 
& Socrates 2003), is likely to be large in the contexts consid- 
ered. As in the work of e.g., Dzhalilov et al. (1992), Zhugzhda 
et al. (1993), and Bogdan et al. (1996) yet more precise de- 
scriptions of the transport should be explored, as should the 
dependence of the stability properties on the temperature and 
density dependence of the opacity (e.g., Bisnovatyi-Kogan 
& Blinnikov 1979; Zhugzhda et al. 1993; Blaes & Socrates 
2003) and the explicit dependence on the scattering albedo 
(e.g., Kaneko & Morita 2006). 

Gradients in the background state — and, particularly, a 
background flux — have been neglected in the analysis of 
^2] Ledoux (1951) considered a consistent background state 
without invoking the Jeans Swindle in calculating the Jeans 
instability and found only a small quantitative change to the 
stability properties of the medium: the Jeans length was in- 
creased by a factor of V2- A detailed assessment of such 
gradients in the context of radiation pressure dominated self- 
gravitating media, as well as an exploration of magnetic fields 
and the physics of the photon bubble instability (Arons 1992; 
Gammie 1998; Blaes & Socrates 2001, 2003) are saved for a 
future effort. 

3. DISCUSSION 

The analysis of §|2] indicates that under the assumption of 
optically-thick equilibrium radiative diffusion, an isotropic 
self-gravitating medium is unstable if the classical gas-only 
isothermal Jeans criterion is not satisfied. If diffusion is rapid, 
the instability is dynamical — the classical Jeans instability. 
If diffusion is slow, then on scales larger than the adiabatic 
Jeans length the medium is dynamically unstable — again, 
the Jeans instability. However, on scales smaller than the 
adiabatic Jeans length the medium is unstable to a diffusive 
mode that acts on a timescale of order the Kelvin-Helmholz 
timescale (eq.[33ll, longer than the dynamical timescale of the 
medium. 

Depending on whether or not the medium is rapidly or 
slowly diffusing (cf. eqs.[39l&l40t. the growth time for gravi- 
tational instability may be significantly decreased with respect 
to the dynamical timescale corresponding to the average den- 
sity of the medium. Here I discuss several astrophysical envi- 
ronments where this analysis is applicable and where it pro- 
vides some insight into the stability properties of the medium. 

3.1. Radiation in Starbursts, ULIRGs, pc-scale AGN Disks, 
& Extreme Massive Star-Forming Regions 

Starburst galaxies are marked by high radiation energy den- 
sity, high gas density, and optical depths to their own dust- 
reprocessed infrared radiation that exceed unity. Scoville et 
al. (2001), Scoville (2003), and TQM have argued that radi- 
ation pressure may dominate the dynamics in these systems. 
Typical temperatures are in the range T ~ 50-200 K and den- 
sities range from n ~ 500- 10 3 cm -3 (e.g., the central regions 
of M82 and NGC 253) to n > 2 x 10 4 cm" 3 (e.g., the nuclei of 
the ULIRG Arp 220; Downes & Solomon 1998) on -lOOpc 
scales. In extreme massive star-forming regions temperatures 
are similar to those in starburst galaxies and ULIRGs, but the 



typical density of clumps and cores responsible for star for- 
mation can be higher. An example is the core of NGC 5253 
with n sw 10 7 cm" 3 and a physical scale of order 1 pc (e.g., 
Turner, Beck, & Ho 2000). 

Yet more extreme physical conditions are expected to ob- 
tain in the self-gravitating pc-scale disks or obscuring "torii" 
thought to attend the process of fueling active galactic nu- 
clei (e.g., Pier & Krolik 1994; Goodman 2003; Sirko & 
Goodman 2003; TQM). There, one expects central disk tem- 
peratures approaching the sublimation temperature of dust 
grains, T su b ~ 10 3 K. Although the gas density is uncer- 
tain in these environments, we can make an estimate by as- 
suming that the disk is marginally Toomre-stable such that 
n w n 2 /(27rGm p ) sw 6 x 10 8 M 8 /?7 3 pc cm" 3 , where M 8 =M/10 8 
M and R\ pc = R/l pc, within the sphere of influence of the 
central supermassive black hole. 4 

Simple estimates indicate that the dust and gas in these sys- 
tems are collisionally and energetically coupled and that the 
IR optical depth is larger than unity. Ignoring the enhance- 
ment of dust-gas coupling due to magnetic fields and grain 
charging, the mean free path for momentum coupling is 



Xd g ~ 10 3 n 4 ' pc : 



10"V pc, 



(41) 



where n x = n/10 v cm . The medium is optically-thick to the 
dust-reprocessed IR radiation field on scales larger than 



A r= i = {up) 



«2.5«4 P CI 



x 10"% 



pc, (42) 



where K2.5 = ft/2.5 cm 2 g" 1 is a representative Rosseland- 
mean dust opacity for T « 100K, assuming solar metallicity 
and a Galactic dust-to-gas ratio (e.g., Figure 1 from Semenov 
et al. 2003). For temperatures near the dust sublimation tem- 
perature, K2.5 ~ 1 is also a fair order-of-magnitude approx- 
imation (but, see Chang et al. 2006). Because \ T= \/\d g ~ 
8 x 10 3 /k2.5 S> 1, the dust and gas are always highly col- 
lisionally coupled if the average medium is optically-thick. 
A rough estimate of the optical depth in the nuclei of the 
ULIRG Arp 220, where the scale of the system is R ~ 100 
pc, is r IR « 10k2.5«4^ioo P c- ^ n systems like M82 and NGC 
253, n is lower and the medium is only marginally optically- 
thick on lOOpc scales. Although observations indicate that 
the obscuring material surrounding AGN may occupy a large 
fraction of Air, theoretical arguments suggest that most of the 
gas may be confined to a thin disk with vertical scale height 
h < R (e.g., Thompson et al. 2005; Chang et al. 2006; Kro- 
lik 2007). The vertical optical depth in such a disk is then 

T IR W 10 2 K 2 .5"8(V0.1pc). 



Although the dust-gas 
on the scales of interest 
perfect. However, in re; 
diation is optically thick 
differences between the 
dense molecular clouds 
tained by a combination 
cosmic rays, and coolin; 



fluid is highly collisionally coupled 
the energetic coupling may not be 
gions for which the cooling line ra- 
we do not expect large temperature 
gas and dust, similar to the case in 
where the gas temperature is main- 
of heating by dust-gas collisions and 
g in molecular lines (see, e.g., Gorti 



lL 4 g T 3 2 pc, where Tj, 



The sublimation radius for dust is S su (, 
7s U b/10 3 K and L46 = Lbh/10 46 ergs s _1 is the Eddington luminosity for a 
10 8 M black hole. 

5 Here, I have assumed a dust grain density of 3 g cm~ 3 and an average dust 
grain radius of aj Ra 0.1 fim. Because grain charging and magnetic fields are 
likely to be important the dust and gas may be regarded as a single, coupled 
fluid on scales larger than Xj g . 



& Hollenbach 2004). At high density, inelastic dust-gas colli- 
sions likely dominate gas heating. Assuming order-unity dif- 
ferences between the gas and dust temperatures, the gas heat- 
ing timescale is roughly 



thesLt/tdy 



:0.03h 4 1/2 7 2 I/2 : 



io-V /2 r- 1/2 



(43) 



where here t& yn = (Gp)~ l l 2 . Assuming tight dust-gas coupling 
and tir > 1, the medium is highly radiation pressure domi- 
nated: 

p r /p g « 10 3 2iV « 10 2 r 3 V . (44) 



Although these estimates imply that starbursts and AGN 
disks are optically-thick and potentially modestly to very 
strongly radiation pressure dominated, the ratio 6 



X ^ (47rGp)'/ 2 
<& s ~ (4p r n/c) 



1400K2.5" 4 /2 7r 



VAr-4 

2.5 n 8 



UkoWT 3 



(45) 



shows explicitly that diffusion of radiation in starbursts is very 
"rapid." The stability properties of these media on the scale 
of the system are thus best represented by the far right-hand 
portion of the left panels of Figures [T] and gas pressure is 
negligible and radiation pressure is important as measured by 
the adiabatic Jeans number $ s , but x ^ $s ar, d tne medium is 
dynamically (Jeans) unstable. 

For fiducial parameters, the scaling for AGN disks also in- 
dicates that they are rapidly diffusing. However, the disk 
parameters in this regime are quite uncertain. For exam- 
ple, Chang et al. (2006) advocate K2.5 ~ 20 or larger for 
a gas density n ~ 10 8 cm" 3 and solar metallicity, implying 
x/$s ~ 0.7r 3 " 4 . The linear dependence of k on metallicity 
and the strong temperature dependence of x/$s implies that 
if grains persist for T > 10 3 K and/or the composition of the 
disk is super-solar, the medium may transition to slowly dif- 
fusing and the timescale for gravitational instability will be 
increased to ?kh (eqs. l33l & 1371 ). Also, an estimate of fir 
shows that gas pressure becomes important on small scales in 
AGN disks and may stabilize the medium in the linear regime 
(TQM). 

3.1.1. The Non-Linear Outcome 

Consider an initial hypothetical equilibrium configuration 
for a self-gravitating disk with starburst/ULIRG-like charac- 
teristics such that {>t ^ lj *&s ~ 1 on ~ 100pc scales and with 
c 2 ss c 2 r w p r /p ~ (JiQ) 2 , where h is the disk scale-height 
(as in the models of TQM). Because diffusion is "rapid," the 
configuration is dynamically unstable on all scales larger than 
the classical gas Jeans length (<C h) and vertical hydrostatic 
equilibrium cannot be maintained in the linear regime. 7 

The non-linear outcome of the Jeans instability in such a 
system is highly uncertain in part because it is tied to star 
formation, which, in turn, determines the character of the ra- 
diation field. The question of whether or not hydrostatic equi- 
librium can be maintained depends crucially on the non-linear 
coupling of the radiation and the gas. One possibility is that 
the large-scale radiation field produced by star formation is 
coupled to the generation of turbulence, which regulates the 
structure of the galaxy and its stability properties. Super- 
sonic turbulence on large scales has been shown to inhibit 

6 Here, (A/Cv) in \ is tak en as K. 1 (see eqs. (20), (26)). 

7 The equilibrium imagined is likely also unstable to convective, magneto- 
rotational, and photon-bubble instabilities (see Blaes & Socrates 2001, 2003). 
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the Jeans instability and gravitational collapse (e.g., Klessen 
et al. 2000; MacLow & Klessen 2004). It has also recently 
been invoked as a basis for understanding the origin of the 
Schmidt/Kennicutt laws (Kennicutt 1998; Krumholz & Mc- 
Kee 2005). Indeed, the turbulent velocities inferred in lo- 
cal starbursts and ULIRGs are large enough that a "turbu- 
lent Jeans number," -d tm \, ~ 6v 2 k 2 /(4%Gp), analogous to -d s 
and $7-, may indicate marginal stability: i9 tul -b « 1. Thus, if 
radiation pressure forces can generate turbulence, perhaps a 
statistical hydrostatic equilibrium can be maintained. 

There are at least two reasons why — in the absence of en- 
ergetic input to the ISM from stars (e.g., supernovae) — that 
star formation may be coupled to the generation of turbulence. 
First, because in the initial fictitious equilibrium state envi- 
sioned radiation pressure is large enough that p r / p w (hil) 2 , 
the radiation field is capable of driving mass motions with ve- 
locities of order Sv w hil if order-unity spatial variations in 
the radiation field are present. Second, although I have not 
shown it in this paper, one expects the astrophysical environ- 
ments described here to be subject to the self-gravitating ana- 
log of the photon bubble instability, which in its non-linear 
state will drive turbulence (Turner et al. 2005). The latter is 
particularly interesting because it motivates a dynamical cou- 
pling between the turbulent energy density (M t urbX me photon 
energy density (M p h), and the magnetic energy density (ug). 

The second of these connections, between M p h and ub in 
galaxies, can be motivated phenomenologically. Recently 
Thompson et al. (2006) have shown that the magnetic field 
strengths in starbursts significantly exceed estimates derived 
from the "minimum energy argument." In addition, they show 
that ub must be a constant, order unity, multiple of u^ in these 
systems (see Condon et al. 1991). This conclusion follows 
from the linearity of the FIR-radio correlation, the radio spec- 
tral indices of star-forming galaxies at GHz frequencies, and 
the fact that the ratio us/uph measures the importance of syn- 
chrotron versus inverse Compton cooling of the cosmic ray 
electrons and positrons (e.g., Condon 1992). The fact the 
galaxies that comprise the FIR-radio correlation have M p h's 
that span five to six dex, and that in the Galaxy M p h ~ m#, 
implies that ub must increase with M p h, from normal Milky 
Way-like galaxies to ULIRGs. The necessity of this lock-step 
increase in both u p h and Ub may signal a dynamical coupling 
between the radiation field and the magnetic field in galaxies. 
Thus, the fact that M p h, ub, and u tuI b are of the same order of 
magnitude may not be a coincidence, but instead a necessary 
consequence of the dynamical coupling between the radiation 
field and the self-gravitating magnetized ISM. 

3.2. Cosmic Rays in Normal Star-Forming Galaxies & 
Clusters 

Although no attempt is made here to model the diffusion 
of cosmic rays, it is instructive to consider the various param- 
eters governing gravitational stability in the case of cosmic 
rays vis a vis radiation. 

The total pressure in cosmic rays in the Galaxy is /? cr ,MW ~ 
10~ 12 ergs cm" 3 (e.g., Boulares & Cox 1990), comparable to 
the energy density in starlight, magnetic fields, and turbu- 
lence. The cosmic ray lifetime is inferred to be t CI w 2-3 x 10 7 
yr (Garcia-Munoz et al. 1977; Connell 1998). Interpreted as 
a diffusion timescale on kpc scales, one infers a cosmic ray 
scattering mean free path of order / m f p « 0. 1 - 1 pc. Addition- 
ally, from the observed grammage traversed by cosmic rays 



in the Galaxy, one infers an average gas density encountered 
by the cosmic rays of n « 0.2 cm" 3 (Engelmann et al. 1990; 
Webber et al. 2003). Writing l m{p = (up)' 1 (cf. eqs. ~ 
also Kuwabara & Ko 2004), tf s sw 50 Pnn^X^L and 8 



0, 



4/5cr 



(4irGp) 



1/2 



0.4/ .lpc.P cr 1 i MW' l 2 2 - ( 46 ) 



Although the parameters are uncertain, equation j46l ) indi- 
cates that cosmic rays are marginally slowly diffusing in nor- 
mal star-forming galaxies. Thus, as for photons in the dense 
pc-scale AGN disks discussed in jj3.ll perhaps on kpc scales 
x/fts may be somewhat less than unity so that the growth rate 
for the gravitational instability is (cf. [37l i 

flOtAdyn ~ 3p cri Mw'o.lpc w 0.2 • 

Depending on the phase of the ISM considered, the isothermal 
Jeans number #r may be very close to unity so that the crit- 
ical x below which the Jeans instability is suppressed can be 
made significantly larger (eq. l32l ). Indeed, cosmic rays have 
recently been proposed as an important large-scale feedback 
mechanism in star-forming galaxies (Socrates et al. 2006). 

Similar estimates may be written down for the central re- 
gions of galaxy clusters, but it is unclear if these regions 
may plausibly be cosmic ray pressure dominated. If they are 
at least modestly so, scaling from equation (f46b for higher 
pressures and lower densities, we see then that for 0.1 < 



mfp 



< 



100pc they are plausibly in the "slow" diffusion limit 
(x/fls < 1); again, the dynamical Jeans instability is quelled 
by the non-thermal pressure support. In this limit, the diffu- 
sive instability identified in equation d32b still acts on ?kh, but 
this timescale is likely many times the age of the universe: for 
c s , r w 1000km s" 1 and n « 10" 2 , f KH « lO^/n^/pc) yr. 

3.3. Individual Massive Stars 

Individual massive stars are radiation pressure dominated 
and slowly diffusing and may in principle also be subject to 
the diffusive mode identified in $2] This is simply a secular in- 
stability of the kind discussed in, e.g., Hansen (1978), and ref- 
erences therein. If present, the growth timescale is the Kelvin- 
Helmholz time, of order ?kh ~ 3kc 2 r /(4irGc) ~ 10 3 yrs for 
typical parameters, where k is the Thomson opacity and c s , r 
is the adiabatic radiation pressure dominated sound speed of 
the fluid ~ (GM / R) 1 1 1 . Although massive stars are known 
to be globally secularly unstable on the Kelvin-Helmholtz 
timescale, it is possible that otherwise stably-stratified (radia- 
tive) regions of their interiors may be locally unstable to a 
variant of the diffusive instability in equation d33l . 

4. SUMMARY & CONCLUSION 

I consider the physics of gravitational instabilities in the 
presence of dynamically important radiation pressure and ra- 
diative diffusion. I find that the medium is always stable on 
scales smaller than the gas-only isothermal Jeans length, Aj.r 
(eq. l23l ). For scales larger than Aj,r there are two pos- 
sibilities depending on whether the me dium is "slowly" or 
"rapidly" diffusing, as defined in 32.51 When diffusion is 
rapid, radiation leaks out of a perturbation without providing a 

8 For simplicity, here I take A/Cv = 1 in %. This is a rough approximation 
in the context of the Galaxy because, depending on which gas phase of the 
ISM is being considered, i9r may be the same order of magnitude as # s . 



sufficient restoring pressure force and the medium is dynami- 
cally unstable on all scales larger than Aj.r, regardless of the 
dominance of radiation pressure. 

The limit of slow diffusion is more interesting. Here, the 
medium is unstable to a diffusive mode at an intermediate 
range of scales between the gas-only isothermal Jeans length 
Aj, t and the larger (gas + radiation) adiabatic Jeans length Aj. s 
(eq. l25l ). The characteristic growth timescale is longer than 
the dynamical timescale. Neglecting gas pressure, it is given 
approximately by equation ( f33l > (see also eqs. |32l & 1371 ), 
which is simply the Kelvin-Helmholz timescale for a radia- 
tion pressure supported self-gravitating fluid parcel to radiate 
its total thermal energy at the Eddington limit. Note that on 
small spatial scales, the characteristic timescale is indepen- 
dent of scale. For A > Aj s the medium is dynamically "Jeans" 
unstable, as expected. Thus, even when radiation pressure 
is dynamically dominant, on precisely the scales where the 
medium is dynamically stable by the usual Jeans criterion 
(A < Aj (S ) it is unstable to a diffusive instability that oper- 
ates on the Kelvin-Helmholtz time. I conclude that radiation 
cannot formally stabilize a self-gravitating medium on scales 
larger than the gas-only isothermal Jeans length. See also the 
discussion of Kaneko & Morita (2006). 

In 33.11 I consider the importance of the results derived in 
^2] for extreme sites of massive star formation including star- 
burst galaxies and pc-scale AGN disks. I argue that the av- 
erage medium in these systems is likely to be radiation pres- 
sure dominated and optically-thick. Importantly, for fiducial 
parameters the photons in these systems are in the rapidly dif- 
fusing limit (xMs > 1; eq- B31 ). For fairly extreme choices 
for the uncertain physical parameters in pc-scale AGN disks 
(e.g., the opacity k) this environment is marginally slowly dif- 
fusing and thus the stability properties of the medium might 
be qualitatively different from rapidly-diffusing starbursts. 

TQM developed a theory of marginally Toomre-stable star- 
burst and AGN disks supported by feedback from radia- 
tion pressure. Because the IR photons produced by dust- 
reprocessed starlight in these systems diffuse rapidly and be- 
cause their characteristic sizes are much larger than the classi- 
cal gas-only Jeans length, the analysis presented here dictates 
that they cannot be supported in the linear regime by radiation 
pressure alone. One may wonder, then, why the entire mass 
of gas in starbursts does not fragment into stars on a single 
dynamical time, in apparent contradiction with observations 



(e.g., Kennicutt 1998). If radiation pressure is to be the domi- 
nant feedback mechanism, then the answer must be that these 
forces are coupled to the generation of supersonic turbulence, 
which may mitigate against complete collapse and fragmenta- 
tion on scales larger than the gas-only Jeans length. In 33.1.11 
I argue that the generation of turbulence likely proceeds from 
the non-linear coupling of the Jeans instability with the ra- 
diation field through star formation, and may be driven by 
the self-gravitating analog of the photon bubble instability. 
This may help explain the apparent order-of-magnitude equiv- 
alence between the radiation, magnetic, and turbulent energy 
densities in starburst systems. Thus, it is important to empha- 
size that the conclusion that radiation pressure alone cannot 
stave off gravitational instability in the rapidly diffusing limit 
does not necessarily imply that the disk cannot be maintained 
in global hydrostatic equilibrium in an average sense by radi- 
ation pressure in the non-linear, turbulent regime. 



In 33.21 I consider the case of cosmic rays diffusing in the 
Galaxy and the cosmic ray halo, and in galaxy clusters. Al- 
though no attempt is made here to calculate the physics of 
cosmic ray diffusion and their thermal coupling to the gas, 
they provide a useful point of contrast with radiation because 
of their very high scattering optical depths. Even so, this con- 
stituent of the ISM of the Galaxy is at the border of "slow" 
and "rapid" diffusion (x/"& s ~ 1) outlined in $2](eq. |46l ) for 
fiducial parameters. In the cluster context diffusion is likely 
more fully in the "slow" limit, but it is unclear if cosmic rays 
dominate the total pressure budget in the central regions (e.g., 
Guo & Oh 2008). Finally, I also briefly mention the possi- 
bility that individual massive stars may be locally unstable to 
this diffusive mode on the local Kelvin-Helmholz timescale in 
otherwise stably-stratified radiative regions of their interiors. 
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APPENDIX 

UNIFORM ROTATION 

Chandrasekhar (1954) and (1961) explored the effect of uniform rotation on the Jeans instability and found that the dispersion 
relation is modified by the Coriolis force in the rotating frame. In particular, he showed that for the special case of waves 
propagating at right angles to the direction ft (cos 2 6 = 0) that equation (Q3 becomes u> 2 = 4H, 2 +c 2 k 2 -4TrGp — that is, if O 2 > TrGp, 
then the Jeans instability is stabilized for any gas sound speed c g (Sp g = c 2 5p assumed). For the general case cos 2 # /0 — for 
waves whose wave vectors have arbitrary angles with respect to the spin axis — if the classical gas Jeans criterion (eq. Q) 
indicates instability (-dj < 1), then the medium is unstable for any |0|. In order to gain some intuition and to make contact with 
the work of Chandrasekhar (1954) it is useful to consider the Jeans instability including uniform rotation, radiation pressure, and 
radiative diffusion. The equation expressing conservation of momentum in the rotating frame is 



+v .Vv = --V/ , -Vi>+2(vxfi). 

at p 



(Al) 



All other equations in the original analysis of §|2]are unchanged. I take k = (0,0, k : ) and fl = (0,Oy,O z ). The perturbation 
equations in component form are 

— ito5p + ipk z 8v z = 



10 

— iu>5v x — 26v y fl z + 25v z Q y = 

—ituSvy + 2Sv x fl z = 

-iu)Sv z + ik z 8P/p+ ik z S$ - 25v x Q y = 

-iujSU + (U+ P)ik z 8v z + ik,5F z = 

-k 2 S<S>-4nG5p = 

ik z 8u r + 3np8F z /c = 0. 

(A2) 

The resulting dispersion relation is (cf. eq. Q3D) 

uj 5 + iluuj 4 - u? [AVt 2 + c 2 k 2 - 4ttG P ] - iujLu 2 [Ati 2 + c 2 T k 2 - 4nGp] 

+ u>(4n 2 cos 2 8) [c 2 k 2 -4TrGp] + iCJ(4Q 2 cos 2 9) \c\k 2 -4-nGp\ =0, (A3) 

where cos 9 = S\/|f2| and Co is the radiation diffusion rate given in equation (fTTT i. Defining 

Q = n 2 /(7rGp), (A4) 

and using the definitions for dj, $s> X> an d €> i n equations (f22b — (|26l>. equation ( |A3t can be rewritten as 

e 5 + ixt? - e 3 (e+^-i) - ix?(Q+$T-i) + eecos 2 0(^-i) + / x ecos 2 0(^ r -i)=o, (A5) 

For Q = 0, equation dA51 > reduces to equation d27l ). In addition, for the special case cos 2 9 = equation (fA5) reduces to equation 
d27| i with the substitutions §r — * Q + "&T and $ s — » Q + "d s . Therefore, for Q > 1 and cos 2 6> = 0, the medium is stabilized for any 
X, $t, and #,, As in Chandrasekhar (1961), I find that for all cos 2 9 =/ 0, if dj < 1, the medium is unstable. 

MORE GENERAL TREATMENTS OF RADIATION TRANSPORT 

The prescription for radiation transport in §|2] makes several approximations. In particular, it neglects the time-dependence 
of the radiation field and it assumes perfect radiative equilibrium so that the radiation and gas temperatures are identical. The 
latter assumption is particularly suspect when diffusion is rapid on a scale k , since radiative equilibrium may not be possible to 
maintain. In fact, contrary to the results of §|2] when the radiation and gas temperatures are distinguished, the gas acoustic speed 
in the limit of rapid diffusion should be the adiabatic gas sound speed and not the isothermal gas sound speed (e.g., Mihlalas 
& Mihalas 1984). More detailed treatments of radiating flows without self-gravity may be found in Dzhalilov et al. (1992), 
Zhugzhda et al. (1993), and Bogdan et al. (1996). Kaneko & Morita (2006) provide a detailed treatment of the radiation that 
distinguishes between scattering and pure absorptive opacity. 

For completeness, here I present an analysis similar to $2] but including the dynamics of the radiation field and allowing for 
energetic decoupling between the radiation and gas. The set of equations is (cf. eqs. (3)-||7)) 

Dp 

Dt F 
D\ 1 

Dt p k ' 

Du„ a 

— - + ju g V -\-Kpc(u r -aT ) = 

Du r 4 , 

— - + -« r V ■ v+ V • F+ Kpc(u r -aT 4 ) = 
Dt 3 

V 2 $-4ttG/9=0 

IDF c 

—fz- + ^u r +K P ¥ = Q, (Bl) 

c Dt 3 

where D/Dt = d/dt + y ■ Vv is the Lagrangian derivative, T is the gas temperature, a = 4<tsb/c is the radiation energy density 
constant, and 7 is the adiabatic index of the gas. The time-derivative of the flux in the Euler equation has been neglected. The 
perturbation equations are 

-iujSp+ipk-6\ = 

-iu>8\ + ik(Sp g /p) + /k<5$ - kSF/c= 

— icuSug + jUgik-8\+ npc(8u r — A g 5T) = 

—iojdu r + (4/3)u r ik-S\ + ik-8F - Kpc(Su r -A g ST) = 

-k 2 S<P-4nG5p = 
-iujSF + (c 2 /3)ikSu r + (npc)5F = 0. (B2) 
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where I have taken 8(aT 4 ) = 4aT 3 ST =A g 5T and employed the Jeans Swindle. The perturbations to the gas pressure and energy 
density are written as 8p g = c 2 5p+(dp^/dT)\ p ST and Su g = c^/Qy— l)Sp + (du g /dT)\ p ST. Solving equations dB21 i. the resulting 
dispersion relation can be written in a number of ways. The form most conducive to comparison with equation dT9b is perhaps 



u 5 + iru 4 (2 + s)- 



r-T.1 1 , i , 2 , 2 ; 2 , ^ u r ,2 

I (1 +s) + v ■ +jcjk H k - 

9 p 



-4irGp 



-/rV 



+r 2 w 



j/*+7Cjr — + 2 + 



4.? 



30 



4M -/t 2 (37-2 + ^)-(2 + s)47rGp 



„2 ,2 



4s 



+ /r^ 2 (4/fc 2 -47rGp)=0, 



+ J — ** (3T-3 + j)-4ttGp I I -•• + -- 

9 p \ 1. - 



where 



4flF 4 



du g /dT\,_ 



v 2 = 



2 1.2 



C l k 



r = Kpc, 



and 



r 2 3 (k^- 1 ) 



K2 



1 1 

3T 2 ' 



(B3) 



(B4) 



The latter is the inverse of the optical depth squared across a scale ~ k , With G = 0, equation dB3b is identical to eq. [101.62] 
of Mihalas & Mihalas (1984). 9 For large optical depth in a radiation pressure dominated medium, ^ 2 /T 2 <C 1 <C s in the third 
and fifth terms of equation ( IB3b - Additionally, for a non-relativistic medium v 2 ^> (4u r /9p)k 2 , c\k 2 , and 4-irGp in the fourth term 
in equation dB3l l. Using this ordering, dividing through by the quantity -T 2 s 2 , and noting that ^ 2 /T = ck 2 /(3np) (cf. eqs. fTTl & 
12610 . equation dB3b becomes 



r u 5 " 


- 


[a; 4 ] 
JTs 


[r 2 * 2 . 



2 . ck j 

+ uj +i- — ur- 



l c 2^ + i^2„ 

3 T 



,.. -47tG P -i— - (c 2 k 2 -4nGp)^0. (B5) 

3k/3 \3 9 p ) inp 

This expression should be compared with equation ( fT9l . Note that the last four terms in equation ( IB5b are qualitatively identical 
to the terms in equation ( fT9] l in the radiation pressure dominated limit (A/Cy) ~ 1. The importance of the first and second 
terms in equation dB5l > are measured by the characteristic frequency Ts with respect to the wave frequencies the expression 
admits. Generically, for sufficiently large Ts/u, these terms are sub-dominant. Thus, if these limits hold, the qualitative stability 
properties outlined in §|2]obtain. 

It is simplest to understand the origin of the extra terms in equation dB3b with respect to equation $1% by taking a step back. 
Neglecting the time-dependence of the flux in the last expression in equation ( IB II ). but leaving the rest of the above analysis 
unchanged, I find that 



u 4 + iTu 3 



-iTu 



, v 



2 , 2 ; 2 , 4 u r , J2 
V S + ^fCjk A IC- 

9 p 



■4irGp 



1,2 



l4k 



As v 2 

— + l + ^r 

3 7 r 2 



4 m,. 



k 2 (3j-3 + s)-4irGp[ l+s + 



r 2 



+ sv 2 (c 2 T k 2 -4nGp)=0. 



(B6) 



The somewhat peculiar terms multiplying ^c\k 2 and (4/9)u r k 2 /p in the fourth term of the dispersion relation (e.g., "4i/37") 
are made clear by examining the explicit and full expression for the sound speed of the radiation and the gas, at constant total 
entropy (cf. eq. |[T8l ). When u r /p ^> c 2 -, one finds that c 2 w (4/3)c 2 ~ + (4/9)u r /p, whereas when u r /p <C c\, one finds that 
c 2 sa 7Cj-+(4/3)(7— l)u r /p. This shows that for large and small s, these terms reduce to the adiabatic sound speed for the radiation 
and the gas, respectively. As in deriving equation dB5t , if I take ^ 2 /T 2 < lc* (second and fourth terms) and v 2 ^> (4u r /9p)k 2 , 

r.2,2 



Cjk , and 47rGp (third term), I find that equation dB6b can be written simply as 



iTs 



rk 2 

3 ■ K 2 
+ U +1- U - 



3 up 



-c 2 T k 2 +-—k 



4u r 2 
97 



ck 2 



4nGp )-i^— (c 2 T k 2 -4TrGp)=0. 
J up 



(B7) 



Compare with equation (lB5b . The importance of Ts is again evident. As emphasized by Bogan et al. (1996) and Blaes & Socrates 
(2003), the characteristic frequency 

[4aT A \ 
Wth = Ts = npc (B8) 

V «« / 

measures the rate at which energy is exchanged between the matter and the radiation field. When this frequency is large, the 
energetic coupling is tight and the analysis presented in $2] is recovered. Thus, for very large Ts and vanishingly small diffusion 



v 2 /r = ck 2 /(3kp) — > — only the second and third terms in equation dB7b survive: 



-u(4c 2 T k 2 



/3 + 



rate across a scale k 

(4/9)u r k 2 / p-4irGp) sw 0. That is, in the optically-thick limit with slow diffusion the radiation pressure contributes to the stability 

of the system against gravitational collapse. As in ^2) the fact that it appears that the system is stabilized if (4/9)u r k 2 /p > 4irGp, 



' Correcting for a sign error in the first term of their eq. [101.58]. 
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even when (4/3)cjk 2 < 4irGp is an artifact of taking the limit of zero diffusion rate. Taking just the last two terms in equation 
dB6l > (the small u> limit) and then taking v 2 jY 2 — ► (high optical depth on a scale k~ l ) and then s — > oo (for tight energetic 
coupling between the radiation and the matter, radiation pressure-dominated), I find that 



u,*-,- 



c z T k l -4nGp 



(4/3)c 2 T k 2 + (4/9)(u r /p)k 2 -4ttG P 



_.ck 2 / 1-tfy 

~ l 3np V(4/3)i? r +i?r-l, '^ 



where 



4 M,fc 2 1 
9 p (4irGp) 

is defined in analogy with i)t and i?,, (eqs. l22l & 12411 ), This expression should be compared with equation ( f32b ; in the radiation 
pressure dominated limit they are identical. Thus, as in {21 1 find that even for highly radiation pressure dominated media, with 
very large optical depth and strong energetic coupling between matter and radiation, the medium is unstable if c\k 2 < 4nGp (the 
isothermal Jeans number d T < 1). As before, at high-fe, the characteristic timescale for instability is independent of spatial scale 
and is simply the Kelvin-Helmholz timescale (cf. eqs. 1331 & 11371 ). 

Alternatively, taking the limit Y — ► in equation ( |B6l l, so that v 2 /T 2 ^> s, v 2 /T 2 ;§> 1, v 2 /T 2 3> v 2 s equation dB6t becomes 

u?-uj{ 1 c 2 T k 2 -4iTGp)^Q. (Bll) 

Contrary to the discussion of $2] which showed that in the limit of rapid diffusion the gas acoustic mode speed is the isothermal 
sound speed cr, here I find that when the energetic coupling between the radiation and the gas is weak, acoustic modes propagate 
at the adiabatic sound speed 7 1 ' 2 cr, as expected (e.g., Mihalas & Miha las 19 84). This effect was not accounted for in the 
analysis of ^because perfect energetic coupling was assumed. Equation (IB lit shows that in the limit Y — ► 0, the classical gas 
Jeans criterion is obtained and that on scales larger than the Jeans length, the medium is unstable, even if (4/9)u r k 2 / p ^> 4-KGp. 
However, equation (1B1U is somewhat deceiving as it may imply to the reader that the medium can be stabilized in the special 
case Cjk 2 < 4irGp, but jCjk 2 > 4irGp in the Y — > limit. This is false, and an artifact of having taken Y = in obtaining equation 
( 1B1U . Expanding instead to first order in Y, I find the unstable mode is 

^-iYs( C \ k2 - 4 ^)=iYs(^L-), (B12) 

\ 1 c 2 k 2 -4irGpJ \jdT-lJ 

which shows that in the special case c\k 2 < 4irGp, but ^c\k 2 > 4irGp, the medium is unstable. 
Expanding equation dB6l ) in the high-£ limit I find that 



1/2 3r 

us ss ±7 ' crk—i — 
2 



4 u r \ s f j— 1 
9p?J + 3 V~ 



(B13) 



in agreement with Blaes & Socrates (2003) (their eq. [57]). In a non-relativistic medium, the second term in square brackets 
dominates so that ui w ±^f l l 2 cjk-iYs(^f- \)/2~f. Thus, in the high-A; limit gas ac oustic waves are damped by emission and 
absorption, again with characteristic damping rate ~ Ys. Although in equation ( IB 13b I obtain a wave speed equal to the adiabatic 
gas sound speed, as in equation dBlU , the high-fc limit is not identical to the limit Y — ► because the gravitational term, which 
dictates stability /instability, disappears at high k. To see this, I write uj = ± \J~fc\k 2 — 4irGp+ iq in equation ( lB61 l. take only linear 
terms in q, and then expand to first order as Y — * 0. I find that 



2;,2 



4 u r \ s (j-l)cjk 



- ±( ^ 2 - 4 ^ )1/2 -'f [{i£) %w-ZG P \ ■ (B14) 

For large k, equation (IB14b reduces to equation dB 1 3I >. However, on scales where gravity is important, equation ( IB 141 > shows 
that if ^c\k 2 > 4nGp then the damping rate of gravity-modified adiabatic gas acoustic waves is altered from the prediction of 
equation (IB13b . More importantly, we see explicitly that if ^c\k 2 < 4irGp the acoustic mode is unstable. 

Taken together, equations ( IB 121 ) and (IB141 > show that in the limit of small Ys, (1) if both $t < 1 and j-dj < 1, then the medium 
is unstable, and (2) that if i3t < I and ~fdj > 1, then the medium is also unstable. In addition, equation ( IB9b shows that when 
Ys is very large and the energetic coupling between matter and radiation is tight, the medium is also unstable for $7- < 1, even 
if § r ;§> 1 , as in $2] Thus, one expects that the medium is only globally stable on a scale k~ l if the classical Jeans criterion is 
satisfied and the isothermal Jeans number is larger than unity, $y > 1 , 

Full Solution to the Dispersion Relation 
In analogy with $7-, "d s , x> an d "&r (eqs. l22l - ||26l & HB10I ). I define the quantities 

a= = , and (= ttt, (B15) 

3 (4nGp) (4irGp)' ^ {4irGp) 1 / 2 

where ( measures the rate of thermal coupling between radiation and gas in units of the dynam ical timescale. With these 
definitions, the approximate solutions and limits of the previous sub-section are illustrated in Figure IbTI which presents the full 
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solution to equation (IB6l l over a broad range of v/T oc r~j/ 1 for the parameters a = 10 7 , 7 = 5/3, i3 r = 3, $7- = 1/5 (fop left panel), 
■dj = 1/2 (fop right panel), dj = 4/5 {bottom left panel), and $7- = 3/2 {bottom right panel) at fixed scale fc . In each panel, I 
take s = 9(7- l)i? r /i?r so that s = 180, 36, 22.5, and 12 from left to right, top to bottom, respectively. As in Figures Q] and [2] open 
and filled circles show the real and imaginary components of £ = ^/{4-irGp) 1 /2 . Positive complex components indicate unstable 
modes. Note that in each panel the damping rates £ sa — r^ and £ sa — j£ are off-scale at intermediate values of i//I\ 

In each panel, because d r — and, by extension, -d s — is larger than unity, in the limit of large Tfc (large T, small v/T) solutions 
qualitatively identical to those obtained in Figures Q] and [2] are recovered. Thus, the left-hand portions of all panels are similar 
to Figure Q] The qualitatively new feature of this figure is the small-Tfc regions in each panel. In the top two panels both c\k 2 
and jcjk 2 are less than 4irGp {'dj < 1 and j-dr < 1) so that the medium is unstable for any v/T. For sufficiently small v/T, the 
growth rate is sub-dynamical and equal to the inverse of the Kelvin-Helmholz timescale (cf. eqs. J33~l & 1371 ). Here, the diffusive 
instability operates. However, for intermediate values of the v/T, the growth rate is dynamical (£ rs ±/(l - ^r) 1 /2 ), whereas for 
very large v/T, £ « ±/(l - 7$7-) 1 /2 . In these cases, the medium is classically Jeans unstable. 



The bottom left panel is different. Here, $7- < 1, but 7^7- > 1 (see eq. I1B12I ). In this special case, two sets of gravity- and 
diffusion- modified acoustic waves exist: (1) at small v/T, the adiabatic radiation pressure dominated acoustic waves are evident 
(as in the top two panels) and (2) at large v/T, adiabatic gas (only) acoustic waves are also present. Note that in this regime (small 
thermal coupling, small 7>) the medium is still unstable, but the growth rate for instability is rs C(l — $t)('Y&t- 1) _1 (eq. IIB12I ) 
in the special case where both the numerator and the denominator are positive. 

The bottom right panel shows a case analogous to the right panels of Figures [T] and [2] with #7- > 1. Here, the Jeans instability 
is stabilized at all Tfc. At intermediate v/T, £ w ±($7- - 1) I//2 , whereas for large values of v/T, £ sj ±(7^7 - 1) 1 /2 . As in the other 
panels, in the limit of strong thermal coupling between the radiation and the gas (large 7>, small v/T), £ rs ±(i?j- 1) 1/2 . 

Scalings & Applications 

The qualitative differences at small optical depth or poor thermal coupling between the radiation and the matter in the right 
hand portion of each of the panels in Figure |BT] with respect to Figures Q] and [2] do not change the conclusions about most of the 
astrophysical applications discussed in §|3]because fiducial estimates for ( and a are very large: 

S {4ttG P ) 1 / 2 {4irGp) 1 / 2 \ u g ) " 4 2 2 ' 5 8 3 

cVpy-T) }p^\ %] 1/2 (p.} 

l m f p {AnGpy/ 2 \pJ 01pc \pj 

a = _4_ = C ^ll « ^_ „ 7 x 10 6 A2 2 »4' w 7 x 10 8 AI 2 V 

7xl0 8 Ak 2 « _I . (B17) 



and 



4irGp 4irGp GpX 2 

Additionally, the ratio v/T is 

v _ 1 f k\_ 1 ( 2tt \ 

f ~ 3V 2 W / ~ 3 172 V^M/ 



OJ^.sn^Aj 1 « 0.03k 2 .5«8 1a -i 



= ^72 (^) «4xlO-\i P cAi. (B18) 

The first line of each equation shows the scalings for starbursts and AGN disks, while the second line shows the scaling appropriate 
for cosmic rays (cf. eq. 1461 ). For all cases considered, £ is very large and v/T is less than unity. The stability properties of these 
media are thus best represented by the left-hand portion of each of the panels in Figure lBTt the radiation and the matter are tightly 
energetically coupled. 
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FIG. B 1 . — Solution to equation 4B6t for a = 10 7 , 7 = 5/3, I?, =3, <&t = 1/10 (top left panel) and for i? r = 1/2 (top right panel), ■&? = 4/5 (bottom left panel), 
$T = 3/2 (bottom right panel), for a very wide range of u/T = t~^ /3 1 ' 2 at fixed scale lC l . Open and filled circles show the real and imaginary parts of the roots 
£, respectively. Positive complex roots indicate instability. Components of the dispersion relation are labeled for clarity. The left portion of each panel (large tj , 
strong thermal coupling) is qualitatively similar to Figures[T]and]2] In each panel $ r > 1 so that j? s > 1. The bottom left panel shows the special case 'dj < 1 and 
y&T > 1. At very large u/T, the gravity- and radiation-modified adiabatic gas only acoustic wave exists and in this region the growth rate for the Jeans instability 
is suppressed. 
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